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LOWER  BOUNDS  TOR  EIGENVALUES* 
by 

E.  Bright  Wilson,  Jr. 
Harvard  University 
Mallinckrodt  Chemical  Laboratories 
Cambridge,  Massachusetts 


ABSTRACT 

Various  methods  for  the  calculation  of  lower  hounds  for  eigen¬ 
values  are  examined ^including  those  of  Weinstein,  Temple,  Baxley  and 
F’ox,  Gay,  and  Miller.  It  is  snown  how  ail  of  these  can  be  derived  in 
a  uni  fied  manner  by  the  projection  technique.  The  all  ornate  forms 
obtained  for  the  Gay  formula  show  how  a  considerably  improved  method 
can  be  readily  obtained.  Applied  to  the  ground  state  of  the  helium 
atom  with  a  simple  screened  hydroponic  trial  function,  this  now  method 
gives  a  lower  bound  closer  to  the  true  energy  than  the  best  upper 
bound  obtained  with  this  form  of  trial  function.  Possible  routes  to 
|  further  improved  methods  are  suggested. 


‘'•This  research  was  made  possible  by  support  extended  Harvard  University 
by  the  Office  of  Naval  Research,  Contract  Nonr-1866  (14).  Any  repro¬ 
duction  in  whole  or  part  is  permitted  for  the  U.S.  Government. 
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A lthour.h  there  exists  a  considerable  literature^  on  the  subject  of 
lower  bounds  to  eigenvalues ,  these  methods  have  been  relatively  little 
used.  The  earlier  methods  all  required  the  calculation  of  H  ,  H  being 
the  Hamiltonian.  This  is  usually  difficult.  In  addition,  the  lower 
bounds  obtained  have  normally  deviated  from  the  true  energy  much  more 
than  docs  the  upper  bound  H. 

The  tremendous  advantage  of  being  able  to  calculate  both  upper  and 

lower  bounds  and  hence  of  endowing  quantum  mechanics  with  a  genuine 

power  of  prediction  has  spurred  the  search  tor  more  practical  and  more 

efficient  methods.  In  the  last  few  years  a  number  of  promising  new 

2 

approaches  have  been  opened  up.  It  is  the  purpose  of  this  paper  to 
show  how  several  of  these  methods  can  be  related  to  one  another,  and,  in 
consequence  to  demonstrate  a  modified  method  of  considerably  improved 
accuracy . 

The  Expansion  Method 

The  older  formulas  can  all  be  derived  by  an  elementary  expansion 

3 

method,  as  shown  by  Crawford  and  Stevenson.  The  trial  function  4  13 
imagined  to  be  expressed  in  terms  of  the  set  of  true  eigenfunctions  t|»n, 
i.e. 

♦  *  EcnV  Ecn2  '  (1> 

One  then  writes 

I  =  /$  (H  -  a)2  $dT  =  H2  -  2aH  +  a2 

=  Ecn2  <En  -  a)2  *  <Eq  -  o)2  (2) 

if  the  constant  a  is  chosen  so  that 

* 

(Ec  -  a)2  *  (Efl  -  a)2  (3) 

the  En  being  the  true  eigenvalues.  The  various  different  methods  result 
if  appropriate  values  are  assigned  to  o.  In  particular  Kato4  has  shown 
by  a  different  argument,  that  when  only  H2,  H  and  a  lower  bound,  E^,  for 


f 


the  first  excited  state  are  known,  then  the  l«:r,t  lower  hound  for  L  is 

o 

,  5 

Temple 1 s 


H  -  H2 

L  >,  EL  =  — - - 

o  '  o  ,  L 


(4) 


1/’  -  li 
x 


which  can  be  derived  above  by  Riving  a  its  maximum  safe  value 


2a  =  Lq  +  tL 
°  t 


(5) 


and  solving*  for  C^. 


It  is  worth  pointing  out  that  Temple's  formula  can  be  modified  to 
apply  to  discrete  excited  states  as  well  as  to  the  ground  state.  Further, 
although  in  practice  it  has  generally  given  a  worse  approximation  to  the 
eigenvalue  than  does  H,  this  is  not  always  true.  In  fact,  in  a  special 
case,  the  exact  eigenvalue  can  be  obtained  from  an  approximate  trial 
function  4.  This  occurs  when  is  exactly  known  and  only  cQ  and  appear 
in  the  expansion  of  $.  However,  it  may  be  empirically  justified  to  get  a 
better  probable  lower  limit  by  averaging  the  Temple  Eq  with  H. 

If  no  information  regarding  Ll  ^3  available,  Weinstein's6  formula, 

—  —  1 
which  uses  H  for  a,  may  be  used^but  ifv  in  fact  H  is  greater  than 

y(po+  the  result  may  not  be  a  true  lower  limit  to  EQ  .  Ordinarily, 

Weinstein's  formula  is  much  les3  efficient  than  Temple's.  It  cannot  be 

better  than  Temple  and  still  be  a  guaranteed  lower  bound.  Even  a  very 

rough  value  for  E^  can  often  give  a  Temple  bound  considerably  better  than 

the  Weinstein  bound.  If  the  Weinstein  lower  bound  is  averaged  with  the 


*A  better  lower  bound  can  of  course  be  found  by  using  the  true  value  of 
E0  or  a  guess  higher  than  Eq  in  the  expression  for  a,  but  although  calcu¬ 
lations  of  this  sort  are  found  in  the  literature,  they  clearly  beg  the 
question.  If  one  knows  a  better  value  of  E^  in  advance,  there  is  no  use 
in  calculating  a  poorer  value  for  this  quantity. 


-II- 


uppcr  bound  ti ,  weighting  the  latter  by  the  maximum  amount  which  permits 
a  guaranteed  lower  bound  to  result,  the  Temple  value  is  obtained. 

Tno  Projection  Me thou 

In  the  atomic  and  molecular  electronic  problem,  the  interelectron 

repulsion  contributes  a  term  to  the  Hamiltonian  which  is  always  positive. 

The  energy  levels  E°  calculated  with  the  omission  of  this  term  are  there - 

n 

fore  always  below  the  true  energies  and  provide  well  known  but  very 
poor  lower  bounds.  The  method  of  intermediate  Hamiltonians^  i3  to  con¬ 
struct  a  series  of  Hamiltonian  operators  whose  eigenvalues  can  be  calcu¬ 
lated  and  which  lie  between  E°  and  E^.  Further,  by  proceeding  to  later 
and  later  members  of  the  series,  the  lower  bounds  can  be  improved  to  any 
extent  desired. 

One  method  of  constructing  intermediate  Hamiltonians  is  the  pro jec- 
7 

tion  method.  Let 

H  =  H°  +  V,  V  >  0  (6) 

and  choose  a  trial  function  <$,  hopefully  approximating  a  true  eigenvalue 
ill  of  H.  Define  a  special  projection  operator  P  so  that  for  any  function 
X  in  the  domain  of  H, 

PX  =  a<t>,  X  =  a*t>  +  h  (7) 

where  the  coefficient  (which  depends  upon  the  function  x)  is  chosen  to 
ensure  the  vanishing  of  the  integral  <nV$>.  Hence,  from  Eq.  (7) 


<$V x>  =  a  <4V<J»  +  <<Vh> 

(8) 

so 

Px  =  a$  =  <$VX>  $ 

(9) 

The  reason  for  this  choice  is  that  then 

t 

<xvx>  =  <pxvpx>  +  «nvn>  5.  <pxvpx>  *  <XVPX>  (io) 

for  any  suitable  x •  This  is  mathematically  sufficient  so  that  the  eigen¬ 
values  E^  of  the  modified  equation 


when  properly  ordered,  satisfy 


E°  .<  E  v<  E 
n  n  K  n 

Further,  a  forma1  solution  for  the  eigenfunction  x  is 

x  =  -<H°  -  EL)"1  VFx  =  -a(H°  -  eV1  V<f> 


<$V$>a  =  <  <pVy> 


when  E  t  E° , 
n 


Further  progress  can  be  made  by  choosing  *  in  terms  of  a  new  func¬ 


tion  such  that 


$  =  V-1  ( H°  -  EL)  K 


Hence 


X  =  -a£ 


and  the  equation  for  u  becomes 

<f,  (H°  -  EL)  V-1  ( H°  -  C!')  a  =  -<E  (H°  -  EL)  £>  a 
which  has  a  non-vanishing  solution  only  if 


where 


<f,  f  (E  )  C>  =  0 

f  (t.L)  =  (H°  -  EL)  V'1  (H°  -  E1’)  ♦  ( H°  -  E,J) 


Provided  that  Eq.  (17)  has  a  real  solution  and  that  V>0,  E  will  be 

a  lower  bound  to  some  eigenvalue  of  H.  In  particular,0  since  E°  4  j-  E 

— —  .  n  n  n 

if  EL  is  below  r.®  it  cannot  be  a  lower  bound  for  E^  and  hence  must  be  a 
lower  bound  for  EQ. 

The  trial  function  S  can  be  in  the  form  of  a  linear  combination  of 

some  finite  set  of  basis  functions  (w  }  and  the  coefficients  can  be  varied 

n 

to  maximize  E^.  With  d£^  =  0,  this  leads  to 

Ec  <u  f  (EL)  w.  >  *  0.  (13) 

n  n  K 

and  nence  to  the  finite  secular  equation 


det  I  <w  f  (E  )  w,  >  I  =  0 
n  K 


The  typical  clement  can  be  written  in  the  expanded  form,  from  Eq.  (18) 


(V_I)nk  (EL)2  -  (H°V-i  +  V_IH0  t  I)nk  C1,  ♦  (H°V"1H°  *  H°)nk  =  0  (n  ) 


“here  ,  ,,-lv  _  „  ,.-1  ,,  .  .  „ 

(V  )  .  s  <u_  v  vl>  etc. 
dk  n  k 

Another  form  can  be  obtained  from  Eq.  (13)  by  the  substitution  = 

# 

V_1r,  or  (  =  (H°  -  EL)_i  c,  in  Eq.  (16).  This  yields 

<r.  |  V"1  ♦  (H°  -  EL)“1  I  ?>  »  0  (22) 

If  t,  is  expressed  as  a  finite  linear  combination  of  eigenfunctions  for  H° , 
* i.e.  as  . 


C  =  1  a 

s  a 

(23) 

0 

,  9 

secular  equation  is 

*  6st  <1;S  -  eL>"1  1  *  0 

(24) 

"  S5t  (E°S  -  eL)  1  *  0 

(25) 

The  form  for  f(L  )  given  in  Eq.  (16)  can  be  usefully  rewritten  by  re¬ 


placing  H°  by  H-V,  so  that 

f(EL)  =  (H  -  EL)  V*1  (H  -  EL)  -  (H  -  LL)  (26) 

This  form  shows  that  if  C  =  ♦  ,  an  exact  eigenfunction  of  H,  then  E**  =  E 

nT  n 

is  a  solution,  because  of  the  factor  H-E^.  It  further  leads  to  the  form, 
for  normalized  E, 

EL  =  H  -  /[( H  -  Eb)  a2  V*1  dt  (27) 

Clearly  V  in  not  in  principle  limited  to  the  interelectron  repulsion.  The 
method  of  proof  required  that  V  be  positive  definite  and  that  E®  of  H°  = 

H-V  be  greater  than  E^.  Otherwise  V  *  is  quite  unrestricted.  In  practice 
the  choice  is  limited  by  the  problem  of  evaluating  the  integral  and  the 
requirement  of  determining  that  E®  >  E^but  it  is  possible  to  consider 
desirable  qualities  for  V.  Thus,  for  maximum  E*1,  V-1  should  be  small;  i.e. 
V  large,  where  C(H-E)C]  is  large.  At  the  same  time,  V  should  be  small 
where  ^  ,the  first  excited  state  for  H-V,  is  large as  to  keep  E®  >  E*\ 


The  simplest  choice  for  V  is  a  constant  A.  Tin-  maximum  value  for 
A  is  set  by 


or 


E°  =  Ex  -  A  *  E“ 

A  =  eJ  -  Co 


(28) 

(29) 


where  is  a  lower  limit  for  the  first  excited  state  of  H.  This  choice 


leads  to 


f(EL)  =  (H  -  Eq)  (H  -  E^) 


or 


r-  E  _  /  i-  E  JJ 


(fij  H  -  H2)  /  ( E j* 


ID 


(30) 


(31) 


which  is  Temple's"  formula1*,  Eq.  (14).  Hence  all  the  older  formulas  can 
be  derived  from  the  projection  technique  for  intermediate  Hamiltonians. 

Improved  Formulas 

Eq.  (2.1)  has  been  applied0  to  the  ground  state  of  helium  with 

various  trial  functions  of  the  Hylleruas  type  and  was  found  to  give 

bettor  results  than  Temple's  formula.  However,  when  it  is  realized  via 
*  -1 

Eq.  (27)  that  V  is  merely  a  weight  factor,  a  considerable  further  im¬ 
provement  is  easily  obtained. 

One  way  is  to  replace  V  a  l/r12  by  C/r12  where  C  is  a  suitably  chosen 
constant  greater  than  one.  For  a  sufficiently  good  trial  function  +  at 
least,  a  higher  (better)  lower  bound  will  be  obtained  by  increasing  C. 
However,  C  cannot  be  made  too  large  or  the  result  will  no  longer  be  a 
lower  bound.  A  sufficient  condition  is  that  C  be  kept  small  enough  so 
that  the  energy  L*  of  the  first  excited  singlet  state  of  the  operator 


H  =  H  -  C/rl2  =  H°  -  (C-l)/ru 


(32) 


is  greater  than  E^  calculated  with  the  given  C.  Since  C/rJ2  is  a  posi¬ 


tive  operator,  if  the  calculated  Eu  were  really  a  lower  limit,  not  to 
the  ground  state,  but  to  the  first  excited  singlet  state  of  helium,  then 


!  9 

il 
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this  E^  would  have  to  be  greater  than  n*‘  i^fined  above.  If  E^  <  £*,  it 

must  be  also  leas  than  EQjas  desired. 

For  helium,  accurate  calculations  are  available  via  perturbation 

S  » 

theory1  for  the  ls2s  XS  state, Nthe  separate  perturbation  terms 
e°  e'  e"  can  be  inserted  in  the  expansion 

E*  =  e°  -  ( C-l)  e’  +  (C-i)2  c"  -  (C-l)3  e,M  (33) 

to  calculate  a  very  pood  number  for  E^  as  a  function  of  C.  Thus  for 
E*  =  -2.903,  C  s  1.73,  a  safe  value.  Then  Eq.  (26)  (or  its  equivalent 
Eq.  (27))  is  solved  for  E^  using  V  =  C/r  .  With  the  simple  screened 
hydrogenlike  trial  function 


,,  -Z'(r  +  r  ) 

<t>  =  N  e  1  2 


(34) 


the  value  Eu  =-2.947  atomic  units  is  obtained  which  is  below  the  presum¬ 


ably  correct  energy  (-2.9037)  by  0.043  units.  This  can  be  compared  with 
the  error  .056  in  the  best  upper  bound  for  this  form  of  function  (which 
is  obtained  with  a  different  value  of  Z').  Presumably  the  lower  bounds 
with  more  complicated  trial  functions  could  be  considerably  improved  by 
this  means  also,  with  no  greater  labor  than  required  with  V  =  l/r^. 

Still  another  improved  form  can  be  obtained  by  using,  instead  of  V, 

V  =  E°+1  -  e£  *  V  (35) 


in  Eq.  (26)  and 


♦  = 


l 

Z 

n=0 


c 

n 


(36) 


Here  E°  and  refer  to  the  problem  with  Hamiltonian  H-V.  The  secular 
equation, will  have  the  form 

det  |  6..  <E°l+1  -  E?)"1  -  <E°t+1  -  EL  t  V)"1  |  «  0  (37) 

This  equation  was  derived  by  Miller^  in  a  different  way  which  has  the 

advantage  of  releasing  the  restriction  that  V  be  positive  definite.  It 

gave  good  results^  for  excited  states  of  helium. 

&C  4  1.8b  was  used. 

W  -  -  •  <*«*>**•*  _  i  _ 


Conclusion 


The  considerations  above  show  that  steady  progress  has  been  made  in 
developing  lower  bound  methods  of  improving  efficiency.  Further  progress 
can  be  expected.  Thus  the  form  of  V  can  surely  be  chosen  so  as  to  meet 
the  requirements  discussed  above  better  than  does  l/r12. 
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